Meandering instability of curved step edges on growth of a crystalline cone 
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We study the meandering instability during growth of an 
isolated nanostructure, a crystalline cone, consisting of con- 
centric circular steps. The onset of the instability is studied 
analytically within the framework of the standard Burton- 
Cabrera-Frank model, which is applied to describe step flow 
growth in circular geometry. We derive the correction to the 
most unstable wavelength and show that in general it depends 
on the curvature in a complicated way. Only in the asymp- 
totic limit where the curvature approaches zero the results 
are shown to reduce to the rectangular case. The results ob- 
tained here are of importance in estimating growth regimes 
for stable nanostructures against step meandering. 

Keywords: Growth, surface diffusion, models of surface 
kinetics. 



I. INTRODUCTION 

Film growth by Molecular Beam Epitaxy (MBE) is es- 
sentially based on the possibility to control growth on 
a submonolayer level. Usually one aims at surfaces as 
smooth as possible with atomistically sharp structures 
consisting of step edges or nanoscale islands. Growth of 
these structures in MBE is affected not only by stochastic 
but also by deterministic instabilities such as step mean- 
dering and bunching. Since the advent of modern film 
growth techniques these surface instabilities have been 
of theoretical and practical interest 1110] . 

To obtain nanostructures in MBE with desired qual- 
ity one of the fundamental concerns is the stability of 
step edges. Instabilities lead to nonuniform layers so it 
is advantageous to suppress them during growth. In step 
flow growth the basic mechanism and properties of the 
meandering instability were treated by Bales and Zang- 
will (BZ) B on the basis of the classic Burton-Cabrera- 
Frank (BCF) model [Q. The case of circular nanostruc- 
tures has so far received much less attention although it 
is an example for a layered nanostructure |g]. Decay of 
mesoscopic circular stepped structures can be described 
in some special cases of stable step flow with good accu- 
racy using the BCF model f^J^]. Also, recently the decay 
and bunching of a circular crystalline cone in the stable 
growth regime has been studied in detail using a similar 
approach [pi. However, to our knowledge there has been 
no study of the meandering instability in the spirit of BZ 
during growth of circular nanostructures even though its 
importance was noted already a decade ago |g[. It is 
thus of interest to know under which conditions circular 



step edges are stable against meandering. In this report 
we present our study of the morphological instability of 
curved step edges. Our viewpoint is based on the BCF 
model and we generalize the BZ results to the case of 
a circular geometry. We perform a linear stability anal- 
ysis and calculate the corrections to the results for the 
straight steps. As expected, our results reduce to the BZ 
case in the limit of an infinite step radius. 



II. STEP FLOW GROWTH IN CIRCULAR 
GEOMETRY 



We study a model of circular steps which are placed 
concentrically on top of each other. The steps can ab- 
sorb and emit atoms which diffuse on terraces between 
the steps. The terrace j is bounded by the steps at rj 
and rj^i. Assuming the flux of adatoms onto and evap- 
oration from the terraces, the adatom concentration Cj 
on the terrace j obeys the well known form of the BCF 
equation hT 



dt 



-l=D,V'c,-c,/r + F, 



(1) 



where Dg is the diffusion coefficient of an adatom on a 
flat terrace, t is the time scale for evaporation, and F is 
the deposition flux. Assuming that the adatom concen- 
tration relaxes much faster than the step edge moves we 
can assume that the terrace is in a quasi-stationary state 
corresponding to dcj/dt — 0. Mass transport through 
the bulk of the material is ignored. 

The model is now fully specifled with the choice of the 
boundary conditions and the requirement of the mass 
conservation at the step edges. The mass conservation 
implies that the edge velocity is given by M 



Vj - V,, 



or ' or 



R], 



(2) 



where i7 is the atomic area, Vj- , Vj+ are the contributions 
to the step edge velocity due to surface currents from the 
upper and lower terrace, respectively, and Rj is the radius 
of the jth step edge. We assume that the velocities Vj^ 
are related to the deviations of the adatom concentration 
from the equilibrium value M 

(3) 



where cj^ 



Vj± = Vlk± [cj — c^'] 



' exp[rK(i?j)] is the equilibrium adatom 



concentration at the edge, F = flj/{kBT), 7 is the free 



energy/ (unit length), Cq'' is the equiHbrium adatom con- 
centration at the straight step, k{Rj) is the local curva- 
ture of the step in the circular geometry (see e.g. Eq. (7) 
in Ref . |8| ) , and fc_ , kj^. are the attachment coefficients as- 
sociated with the upper and lower terraces, respectively. 
From Eqs. (y|) and (y|) we obtain the mixed boundary 
conditions at the step edge 



Dc ■ 
dc 



^k. 



Cj|fl,-C 



eq] 



Ds-^\B, + ^ =fc-[Cj|fl, + i-C^'%] 



(4) 



Defining a new field Uj = Ci — tF Eq. (|I|) becomes the 
Helmholtz equation in the stationary limit 0,0] : 



2 1 



(5) 



where Xg = j/DsT is the diffusion length. The solu- 
tion of Eq. (|3) for the perfectly circular step is given 
by v3{r) = a°Io(r/xs) -I- b°Ko{r/xs), where Io(x) and 
Kq{x) are the zeroth order modified Bessel functions ||] 
and a° and 6° are coefficients determined by the bound- 
ary conditions. In the linear stability analysis a small 
perturbation is added to the step edge and the equations 
are solved to first order in the perturbation amplitude. 
We set rj{9) = Rj + eexp[m6' -I- ujt] + c.c, where e is a 
small parameter, uj a growth rate, |n| > 1 is an integer 
(due to periodicity), and c.c. denotes the complex conju- 
gate. The solution to Eq. (||) in the first order in e gives 
Uj{r,0) = u°{r) + e[Aj"I„(r/x,) + B^"K„(r/x,)]e™«+'"*, 
where K„(a;) and In (a;) are the modified Bessel functions 
of integer order n, and the coefficients A" and i?" are 
determined by the boundary conditions. When the so- 
lution is found, the growth rate u can be deduced using 



Eq. 



and Vi 



linearly unstable. 






ujfj{9). If w > the step edge is 



III. MORPHOLOGICAL INSTABILITY OF A 
CIRCULAR STEP 

Qualitatively, the growth of stepped structures is un- 
stable against step meandering when the flux of adatoms 
from the upper terrace is reduced, e.g. due to the Ehrlich- 
Schwoebel barrier. This is basically the origin of the mor- 
phological instability on vicinal surfaces (the BZ case) [gj . 
However, in the case of a circular step, the stabilizing 
effect of the step curvature is expected to be more pro- 
nounced than in the rectangular geometry. Therefore, we 
expect the possible instability to be weaker than in the 
rectangular case since the line tension tends to smoothen 
the steps. Here we consider the cases fc+ ^ oo, fc_ = 
(one-sided model), and k^ ^ fc_ non-zero and finite 
(asymmetric model). 



A. One-sided model 

In the one-sided model fc+ ^ cxi corresponds to instan- 
taneous attachment from the lower terrace and /e_ = 
implies an infinite Ehrlich-Schwoebel barrier. In this 
limit the velocity of the step with radius R is given by 
V = V+ = DQ,{du/dr) and the stability function be- 






uj{n) 



hi 



c„-|-(l 
P 



n\ (6) 



where ^s = T/{xstAF) is the capillary length, AF — F— 
Cq'^t, and p = R/xg. The coefficients are given by a„ = 

[iIk„ - I„K^][I;Ko + IqK;], bl = [iIk„ - I„K^][I'iKi - 

fiK'i], hi = [v^w^ - iIk;j[Iik; - i^i^i], hi = [i„k; - 

i;,K„][llK;-i;K'i]^andc„ = [i1k„-I„K^][i],K„-I„K^], 

where I„ = ln{p), In = ^nip+l/xs), I is the terrace width, 
and the prime indicates the derivative with respect to the 
scaled variable p. By using the asymptotic formulae of 
the modified Bessel functions for n, R large with n/R = 
q — const, f^ the growth rate (j^) can be shown to reduce 
to the result of Bales and Zangwill Q in the limit i? -^ oo. 

The stability function is plotted in Fig. 1 and it ap- 
proaches a limiting form when the radius of the step in- 
creases (curvature decreases). As can be seen from the 
figure, there exists a limiting value Re for the radius such 
that steps with R < Re are always stable. For radii 
R > Re there exists a critical value of the wavevector 
Qe such that for q > qe the edge is stable. The criti- 
cal wavevector depends on curvature and with increasing 
curvature qe decreases, i.e. the critical wavelength where 
the instability sets in is shifted to larger wavelengths. 

In the limit of large n, R we obtain the correction term 
to the results of Bales and Zangwill of the order of 1/R 
as 

^^=-Bz(g) + ^E„ (7) 

where ujbz is the rectangular result (Eq. (13) in Ref. [^), 
S, = [xslq^K-^ + A-1 ta.nh{Kqi){Pq^ + A~^) + tanh([) 
-26] X ^ 



cosh(Aq^) cosh(Z) 

p) + (P<z2+A-2)tanh2(y 
+2^s + isX^q^ - ^sXsq\Pq^ + A"') tanh2(A,[) 



- tanh(f)(2 + fq^) + {fq^ + A^^) tanh^(Ag[) tanh(f) 



Aq — -y/l + (x^g)^, and / — Ijxg. Using the equation 
above we can determine analytically the corrections to 
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Fig. 1: The growth rate uj{q) as a function of the 
wavevector q = n/R with a moderate value of the ter- 
race width / in the case of one-sided model. The stable 
regime corresponds to w < 0. The radius of the step R 
increases in the direction of the arrow. The curves in- 
dicate that structures smaller than a critical size R^ are 
stable against step meandering. 
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Fig. 2: The critical wavevector Qc as a function of the 
capillary length <^s in the case of the one-sided model. 
The circles are from of Eq. ffl), the upper solid line is 
the BZ result [pi, and the lower solid line is Eq. (|8|) for 
I ^ Xs- The inset shows the case I <^ Xs, with the lower 
solid line given by Eq. (g). The approximations follow 
the numerically plotted curve, whereas the BZ results 
deviate for larger values of £,s ■ 



the critical wavevector defined as oj{qc) = 0. The results 
are for I '^ Xs'- 



^ sqc 




T + Mt-'^) >2;^9c>1; 



2R \ e 

1(1-26) 



Jj- 1 ^sQc "^ J-! 



(8) 



and for / < x^ (and rq^ <C 1) 



^sQc 




R' 



(9) 



Omitting the \/R terms we obtain the BZ results for the 
rectangular geometry |^. In Fig. 2 the curves for the 
critical wavevector q^ against the capillary length ^ are 
shown ior I ^ Xs, Xsqc ^ 1- The corrected result fol- 
lows the numerically plotted curve whereas the BZ result 
deviates considerably. For ^ large enough the edge is 
always stable and qc aprroaches zero. The inset of Fig. 2 
shows the case I <^ Xs which behaves in a similar way. 



B. Asymmetric model 



QAF 



7j 



7j-i Ij 



P,- 



7j-r 



-]K{Pi) (10) 



where aj , (5j , and 7^ are the coefficients related to unper- 
turbed steps, and A\^^ B-}^, and I?^ are obtained from the 
expansion linear in e. Define d± = d±/ Xs = {Ds/k±)/xs, 

^i,± = ^n{Pj)Td±K'„{pj) andX;^ ± = ln{pj)Td±lniPj)- 
Then the coefficients are given by 

(6/p, - 1)IC^} - (e./p,H 



a 



l)ICl 



rJ+i 



Ps = -iis/p^ - 1)1^ + (^-//'»+l - 1)^0,+ - 






7j = '^i 

-6(n2-l)[/Ci+/p^\i-/C^t' 



/C;-_1+/3,[/Cl 



x- 



i+1 

n,— 



IP 



Bl 



-".^^2^ 



■J + ItJ 



n,+ 



-Ji^:^]-P,llCi,^; 



■J+1 






Hsin" -l)[ll+lp%,-Iit-/p]] 



- _TJ+^ir3 



VI = -Xr_/C: 



T^ AT 



J + 1 



When the kinetic coefficients at the step edge from 
the lower and upper terrace are both finite and non- 
zero, the velocity of the jth step edge is given by Vj = 
Ds^[{duj/dr) — (9Mj_i/9r)]|r=fj(e) and the growth rate 
becomes 



The resulting expression for the growth rate ijj{q) is 
rather complicated and we have not found any essen- 
tially simpler form which would be accurate enough for 
all curvatures. However, the essential qualitatitive fea- 
tures can be extracted from the special cases, and the 
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Fig. 3: The instability function uj{q) as a function of 

the wavevector q in the asymmetric case. The attach- 
ment coefhcients fc_, A:_|_ approach the same value from 
top to bottom. For fc_ < A:+ the steps are always stable. 
The inset shows lu with increasing phase difference (j) be- 
tween two adjacent steps. Increasing (jj makes steps more 
stable against meandering as in the rectangular case pOl . 



behavior is similar to the one-sided case. The numerically 
plotted values of w indicate that as k- approaches fc+ the 
growth rate w < at all values of q. This is shown in 
Fig. 3. We can thus conclude that the one-sided model 
is the most unstable. 

So far, we have only considered the case of in-phase 
step growth. The stability of growth depends also on 
the phase of the neighboring steps, and the above anal- 
ysis can be extended to the more general situation with 
an arbirary phase between adjacent steps. However, the 
general results give only little additional insight and only 
the numerical results are shown in the inset of Fig. 3. In- 
creasing the phase difference between the two adjacent 
steps makes the stegs more stable as in the case of rect- 
angular geometry llfl] . 



terraces to the step are both finite and non-zero behaves 
qualitatively similarly to the one-sided model. 

The results presented here extend the results of the BZ 
model for rectangular geometry and they are of interest in 
modeling e.g. the behaviour of isolated crystalline cones. 
Results given here can be used to check whether the sim- 
plified model which assumes the perfectly circular step 
edges in decaying cones is valid. If the meandering insta- 
bility is expected to be present, it can play a significant 
role in the evolution of nanostructures. 
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IV. DISCUSSION AND CONCLUSIONS 



In this work we have generalized the meandering insta- 
bility to structures made of concentric islands. We have 
shown that in the case of circular cones growth instability 
due to terrace diffusion can arise as in the case of rect- 
angular geometry. However, the instability is suppressed 
by the curvature and structures with smaller sizes than 
a critical size are stable. The critical size depends on 
the microscopic parameters of the system. We present 
here also the asymptotic analytical results for small cur- 
vatures. In the limit R ~* oo the growth rate approaches 
the expression found by Bales and Zangwill 0] . The more 
general case where attachment from the upper and lower 



